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Abstract. Recent optical spectra of post-AGB stars show the presence of 
C2, CN, and CH + originating in the circumstellar shell. We present here 
new, higher resolution spectra which provide constraints on the physical 
parameters and information on the line profiles. An empirical curve of 
growth for the C2 Phillips and CN Red system lines in the spectrum of 
HD 56126 yields b = 0.50±g;|§ km s -1 . CH+ (0,0) emission lines in the 
spectrum of the Red Rectangle have been resolved with a FWHM ~ 8.5± 
0.8 kms" 1 . The circumstellar CN lines of IRAS 08005-2356 are resolved 
into two separate components with a velocity separation of Av = 5.7 ± 2.0 
km s _1 . The line profiles of CN of HD 235858 have not been resolved. 



1. Introduction 

Post-AGB stars are in a transition stage between the Asymptotic Giant 
Branch (AGB) and the planetary nebulae (PN) stage. During the early 
stage of post-AGB evolution the star is obscured by material expelled dur- 
ing the AGB phase (the AGB ejecta). As this ejecta slowly moves away 
from the central star, the optical depth decreases and the star can be de- 
tected in the optical. When the star reaches high enough temperatures, the 
AGB ejecta is ionized and is observable as a planetary nebula. 

We have studied optically bright post-AGB stars (spectral type A to G 
supergiants) which show circumstellar molecular line absorption (C2 and 
CN, or CH + ) or emission (CH + ) in their optical spectra. The radial ve- 
locities and low excitation temperatures of the molecules (Bakker et al. 
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1996a,b) identify them as circumstellar rather than photospheric or inter- 
stellar (see also Hrivnak 1995). The excitation of C2 is generally described 
by rather high temperatures, T ex = 43 — 399 K, whereas the CN excita- 
tion is found to be much lower, 18 — 50 K. The reason for this difference 
is discussed in section 2.2. The observed abundances and excitation of the 
molecules can lead to a determination of the mass-loss rate. 

We have used the C 2 (A 1 ^ - X 1 ^) and CN (A 2 II - X 2 £+) data of 
HD 56126 to study optical depth effects by means of the curve of growth, 
and we present the first results of a survey to observe these optical molecular 
bands at high spectral resolution (R > 120 000). Our primary goal is to 
resolve the line profiles and to determine the Doppler parameter b and 
the chemical (e.g. abundances) and physical (e.g. expansion velocities and 
temperatures) conditions of the circumstellar shell. 

2. Curve of Growth Analysis for HD 56126 

2.1. THE EMPIRICAL CURVE OF GROWTH 

A curve of growth (CoG) has been empirically determined for the 12 C 12 C 
Phillips (A x n u - X X S+) (v' = 1, 2, 3, v" = 0, J" < 24) and the 12 C 14 N Red 
system (A 2 n - X 2 S+) bands (v' = 1,2,3, v" = 0, N" < 3) of HD 56126 
(Figure 1). Our motive for this investigation is to decide whether the as- 
sumption of optically thin lines for the weaker bands (e.g. 3,0) is allowed. 
The equivalent widths are taken from the work of Bakker et al. (1996a,b). 
The range of line oscillator strengths is 6.67 x 10~ 5 to 1.44 x 10 -3 and 
7.19 x 10" 5 to 9.79 x 10~ 4 for C 2 and CN, respectively. 

Because three C 2 Phillips bands have been observed, the CoG has for 
each lower energy level J" up to nine transitions (a P, Q, and R branch 
for each band). This redundancy allows the determination of the CoG. 
All transitions from a given J" level are fitted to the CoG by changing the 
column density of that level N(J"). For the CN Red system the redundancy 
is higher due to spin-doublet splitting of the lower and upper electronic 
state and A-type doubling of the upper electronic state. There are twelve 
transitions from each A^" level: six main and six satellite branches. Three 
bands were used, which gives in total 36 allowed transitions per N" level 
(Fi and F 2 ). 

After all observed J" (or N") levels have been fitted to the CoG (Fig. 1, 
left panels) an optical-depth-corrected absolute rotational diagram (Fig. 1, 
middle panels) gives the absolute population for each J" (or TV") level. Un- 
der the assumption of a Boltzmann distribution, a linear fit to the diagram 
gives the average rotational temperature T TOt . Finally the rotational tem- 
perature from two successive energy levels can be determined as a function 
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Figure 1. Curve of growth analysis for three C2 Phillips and three CN Red system 
bands in the optical spectrum of HD 56126 (diamond: v' = 3, asterisk: v' = 2, plus: 
v' = 1, all for v" = 0). The upper three panels concern C2 and the lower three panels 
concern CN. Left panels: CoG with the best theoretical CoG over-plotted. Middle panels: 
optical-depth-corrected rotational diagram with a fit to the data assuming a Boltzmann 
distribution (dashed line) and a smooth non-Boltzmann distribution (solid line). Right 
panels: the rotational temperature determined from two successive energy levels (the 
dashed line is the average T TOt ). 



of the lower energy level (Fig. 1, right panels). The parameters derived from 
this analysis are given in Table 1. 

2.2. RESULTS AND INTERPRETATION 

The empirical CoG (Fig. 1) shows clearly that the observations cover the 
optically-thin and saturated parts of the CoG. The theoretical CoG has 
a Doppler parameter b = 0.50^o' 2 3 ^ m s_1 an< ^ t ~ 1 is reached at an 
equivalent width of EW{t ~ 1) = 18 mA. Since the strongest lines in the 
(3,0) transitions of C2 and CN have equivalent width of 38 and 34 mA, 
respectively, the optically thin approximation is only valid for the weaker 
lines in these bands. We note that all the lines of the C2 (y' ,v") = (3,0) 
band seem to be offset with respect to the other lines. This might indicate 
that the band oscillator strength used is somewhat too low. 
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TABLE 1. Results from the curve of growth analysis of C2 and CN in 
the optical spectrum of HD 56126 



C 2 A 1 n u -X 1 £+ CNA 2 n-X 2 E+ 



b 


0.50(+0.59, -0.23) 


0.50(±0.59, -0.23) 


km s" 


T r ot average 


336 ± 34 


12 ± 1 


K 


Tot minimum 


189 ± 10 


9 ±10 


K 


Tot maximum 


1623 ± 10 


15 ±10 


K 


log Nobs 


15.38 ±0.10 


15.56 ±0.10 


p cm 


Whelio 


79.35 ±3.12 


77.78 ±1.17 


km s" 



The right panel of Fig. 1 clearly shows that the rotational temperature 
for C2 is not constant. The molecule is therefore not in local thermodynamic 
equilibrium (LTE) and the population distribution over the rotational en- 
ergy levels is non-Boltzmann. Van Dishoeck & Black (1982) have shown 
that interstellar C2 is radiatively pumped. C2 is a homonuclear molecule 
and does not have allowed pure rotational or vibrational transitions and 
can therefore not cool radiatively: T rot > T^ Q . For low J" levels the ro- 
tational temperature reaches the kinetic temperature, while for very high 
J" levels the rotational temperature is expected to reach the color tem- 
perature of the local radiation field. CN on the other hand can effectively 
cool: T ro t < Tk in . Based on the population ratio between the C2 Phillips 
band J" = and J" = 2 levels we find: T kin = 189 ± 10 K. Combining 
this with the derived b yields f m icroturb = 0.34 ± 0.6 km s _1 . The measured 
b = 0.50l~o;23 km s- 1 gives a line profile with FWHM = 2v / rn2 x b = 0.83 
km s" 1 . In order to resolve these lines a spectral resolution of R > 360 000 
is needed. In the presence of macroturbulence the lines are broader and 
can be resolved at a lower spectral resolution. The microturbulence is very 
likely due to a velocity gradient in the line of sight. 

3. Line Profiles 

3.1. CH+ EMISSION LINES OF THE RED RECTANGLE 

The CH + A 1 !! — X X S + (0,0) emission band in the optical spectrum of the 
Red Rectangle (HD 44179) has been observed at a resolution of R ~ 120 000 
using the 2.7-m telescope of McDonald Observatory (Figure 2). These lines 
originate from levels which are 34 000 K above the ground level. We have 
resolved the line profile of the strongest CH + emission lines and find a 
FWHM « 8.5 ± 0.8 km s^ 1 (Table 2). The intensities of the emission lines 
fall below the detection limit for FWFM rs 20.0 ± 1.0 km s^ 1 . The line 
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Figure 2. Left: normalized line profiles of the strongest CH + A 1 !! — X x £ + (0,0) emission 
lines and CO radio emission lines (after Jura et al. 1995) of the Red Rectangle. The thick 
dashed lines at 17.7 km s gives the average velocity, with on both sides (at a velocity 
offset of 10.0 km s _1 ) a line where the intensity of the emission line has fallen below the 
detection limit. The dashed line at approximately —0.7 km s _1 (J" dependent) marks the 
location where 13 CH + is predicted. Middle: A non-detection of 13 CH + after the signals 
of several lines have been co-added on a velocity axis. It seems that the 12 CH + lines are 
a disk-like feature. Right: the Ha profile with the central spike at the system velocity. 



profiles of the R(l), Q(l), and Q(2) lines suggest a central absorption or the 
presence of two emission components at 16.0 and 18.1 km s -1 . The location 
of 13 CH + is at approximately ~0.7 km s _1 (J" dependent). To investigate 
the presence of weak isotopic lines we have added the signals of all lines 
into one single profile (middle panel of Fig. 2). We have a non-detection of 
13 CH + which yields an isotopic ratio of 12 C/ 13 C > 22. This is consistent 
with the progenitor being a carbon star (Smith & Lambert 1990). 

The line profile of CO is interpreted as due to a broad and narrow 
component (Jura et al. 1995, 1996). The broad component has a width 
comparable to that of the CH + lines which might suggest that they are 
formed in the same region (the circumbinary disk). The Ha and CO profiles 
shows a strong central emission at the system velocity. Jura et al. argue 
that the spike is formed in an extended region of ionized gas. 
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TABLE 2. 


Results for the Red Rectangle, HD 235858, and IRAS 08005-2356 




Red Rectangle 


HD 235858 


IRAS 08005-2356 










pw T}ar\ (9 f)\ 
v^i\ jxea yA,u) 




Date 


1995 Dec. 13 


1995 Sep. 17 


1995 Dec. 12 




HJD 


2450064.7428 


2449977.6607 


2450065.8708 




^*,hclio 


27.7 ± 1.4 


-34.1 ± 0.1 


64.5 ± 1.0 


km s 1 


^mol.hellio 


17.7 ±0.1 


-48.2 ± 0.1 


19.4 ±0.2 (A) 


km s _1 








25.2 ±1.0 (B) 




■ufwhm 


8.5 ±0.8 


2.2 ± 0.1 


3.1 ±0.9 (A) 


km s _1 








3.4 ± 2.0 (B) 




dfwfm 


20.0 ± 1.0 






km s _1 


12 C/ 13 C 


> 22 


> 11 


> 11 





3.2. CN ABSORPTION IN HD 235858 AND IRAS 08005-2356 

HD 235858 has been observed at a resolution of R ~ 120 000. Figure 3 
shows a part of the spectrum which contains three different categories of 
molecular features. The broad absorption lines are due to photospheric CN 
(T e ff = 5500 K), the narrow absorption lines are circumstellar CN (T e g 
= 20 K and logiV = 15.30 p cm -2 ), and the three strongest features arc 
telluric H2O. There is one photospheric atomic line (Nl) present in this 
spectrum. The circumstellar CN lines are not resolved. Since Zacs et al. 
(1995) did not notice photospheric CN absorption in their spectrum, it 
seems that this pulsating star has only photospheric molecular absorption 
when the stellar effective temperature is sufficient low (when the star is 
largest). 

In an earlier report on the detection of circumstellar CN in IRAS 08005- 
2356 we reported the presence of only one broad component (Bakker et al. 
1996b). The new spectra (Fig. 3, R « 120 000) have resolved the broad 
component into two separate resolved components with a velocity separa- 
tion of Av = 5.7 ± 2.0 km s -1 . The two absorption components could be 
due to wind material moving at two different velocities (multiple shells), or 
to two photodissociation fronts at different expansion velocities — possibly 
one due to the stellar and one to the interstellar radiation field. 

4. Summary 

We have investigated optical depth effect for the circumstellar C2 and CN 
lines in the optical spectrum of HD 56126. From a curve of growth anal- 
ysis we find a Doppler parameter b = 0.50i~o;23 km s -1 . CN (2,0) lines of 
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Figure 3. Upper: a part of the optical spectrum of HD 235858 which shows the CN 
Red system (2,0) band lines. Lower: for IRAS 08005-2356, a synthetic spectrum using 
photospheric and circumstellar CN lines, and Nl (7898.985 A) is over-plotted (thin line). 
The abscissa is in the rest frame of the molecule. 

IRAS 08005-2356 and HD 235858 and CH+ of the Red Rectangle have been 
observed at a resolution R > 120 000. CH + has been resolved. The CN lines 
of IRAS 08005-2356 are resolved into two separate resolved components, 
and the CN of HD 235858 has not been resolved. 
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